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ABSTRACT 

We report the discovery of two very cold and massive molecular cloud cores in the region ISOSS 
J18364-0221. The object has been identified by a systematic search for very early evolutionary stages 
of high-mass stars using the 170 /im ISOPHOT Serendipity Survey (ISOSS). Submm continuum and 
molecular line measurements reveal two compact cores within this region. The first core has a tem- 
perature of 16.5 K, shows signs of ongoing infall and outflows, has no NIR or MIR counterpart and 
is massive enough (M ~ 75 M ) to form at least one O star with an associated cluster. It is there- 
fore considered a candidate for a genuine high-mass protostar and a high-mass analog to the Class 
objects. The second core has an average gas and dust temperature of only ~ 12 K and a mass of 
M ~ 280 M© . Its temperature and level of turbulence are below the values found for massive cores so 
far and are suggested to represent the initial conditions from which high-mass star formation occurs. 
Subject headings: Stars: formation — ISM: dust, extinction — ISM: kinematics and dynamics — ISM: 

clouds — ISM: individual objects: ISOSS J18364-0221 — Astronomical data bases: 

ISO data base 



> 

o 

t> 

&\ 
O 

in 
o 

a: 

6 

in 



1. INTRODUCTION 

Most of the current knowledge about the formation and 
early evolution of high-mass stars has been inferred from 
observations of luminous inf rared sources: High-mass 
protostellar objects (HMPO) ( Sridharan et al I2002D and 
hot molecular cores fHMCsl i|Kurtz et alJl200C|L both 
among the youngest known manifestations of high-mass 
stars, are dense condensations of warm (T > 30 K) gas 
and dust. Given their high luminosity of Lir, > 10 3 L 
most of these objects were readily detected by the IRAS 
mission and have been identified by using color cri- 
teria on their infrared fluxes (e.g. jFonta ni el; alj 2005: 
Williams et alJl2004l:lMolinari et al.ll2002HHenning et all 
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Much less is known about a potential evolutionary 
stage prior to the HPMO/HMC phase. It has been sug- 
gested that very cold (T ~ 10 K) and massive cores might 
repre sent the initial con ditions from which massive stars 
form l)Evans et al]l2002l) . similar to the pre-stella r cores 
in the low-mass regime l|Ward-ThomDsonl |2002|) . Al- 
though promising candidates for very young protostars 
have been recen tly discovered near more evolved and lu- 
minous objects ijGarav et alJl2004t iForbrich et alJl2004t 
IWu et alJ2005|) . no clear example of such a very cold core 
is known so far. 

Using the 170 /urn ISOPHOT Serendipty Survey 
(ISOSS) lILemke et alJlToM IBogun et al]ll99ffl we have 
performed an unbiased search for such very cold sources 
based on their Fno^, m / '-Fioo/an color. By cross-correlating 
the ISOSS catalog with the IRAS point source catalog 
and molecular line surveys, we have identified more than 
50 massive [M > 100 Mq) candidate sources with a flux 
ratio of Fivo^m/ -FioOjum > 2, implying average dust tem- 
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peratures lower than 18 K ijKrause et al.ll2004]) . 

We have launched multi-wavelength follow-up obser- 
vations to explore the physical conditions in this unique 
sample of cold an d massive star forming regions (e.g. 
iKrause et al.l l2003) . In this letter we report the discov- 
ery of extremely cold and massive cloud cores in the re- 
gion ISOSS J18364-0221 which provide new constraints 
on the physical conditions under which high-mass stars 
form. 

2. OBSERVATIONS 

Fig-fflshows far-infrared observations of ISOSS J18364- 
0221. The object was centrally crossed by two de- 
tector pixels of the C200 camera onboard ISO (beam 
size of 1.5'FWHM) during an ISOSS scan. A 170 /im 
flux Fi7o (UBl = 235 ± 47 Jy has been derived using 
procedures bv IKrause! l|2003l) and iStickel et alJ {2004). 
The compact ISOSS source coincides in position with 
the unresolved IRAS point source 18339-0224 (R.A. 
(2000)=18 h 36 m 24!7, DEC. (2000)=-02°21'49") which 
has previously be en detected in the CO (1-0) transition 
ijYang et al J 120 02). Its CO radial velocity of «lsr ~ 
33km/s correspond s to a kinematical di stance of d ~ 
2.2 kpc according to lBrand fc BlitJ l) 19931) . This is con- 
sistent with the distance derived by the empirical relation 
d = 320V 57 pc, where V is the number of foreground 
stars detected on t he POSS-B image inside a 5' aperture 
towards the cloud l)Herbst k, Sawverlll981|) . 

Submm continuum jiggle maps at 850 /jm and 450 /im 
were obtained with SCUBA {Holland et al.lll999l) at the 
James Clerk Maxwell Telescope ( JCMT) in May 2003 un- 
der good atmospheric transmission conditions (T225GHZ = 
0.067±0.003). The data was reduced with SURF, includ- 
ing sky noise reduction, yielding a photometric accuracy 
of 20% (450 /xm) and 12% (850 /im). The beam size is 
7.5" and 14.5", respectively. 

Observations of the Ammonia inversion lines (J,K) = 
(1,1), (2,2), (3,3), and (4,4) were carried out at the Effcls- 
berg 100 m telescope in January 2003. The beam size at 
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the observing frequency of 23.7 GHz is 40". The pointing 
accuracy is estimated to 5" rms and the flux uncertainty 
to 10%. 

Molecular line measurements of CO(2-l) and C 18 0(2- 
1) (beam size - 11"), HCO+(3-2) and H 13 CO+(3-2) 
(beam size ~ 9.5") were obtained at the IRAM 30 m in 
November 2003 and of CO(3-2) (beam size ~ 22") at the 
Heinrich Hertz Telescope in December 2004. We used 
the facility receivers and reduced all spectra with the 
CLASS software. The atmospheric transmission wwas 
(r 2 25GHz ~ 0.15) for the HCO+(3-2) (r 225GHz ~ 0.5) for 
the CO(2-l) and (t 22 5ghz ~ 0.1) for the CO(3-2) obser- 
vations. The estimated pointing error for both telescopes 
is 3" and the radiometric accuracy 15%. 

Near infrared images in J, H, and Ks have been ob- 
tained at the Calar Alto 3.5 m telescope in Oktober 2003 
and June 2004 with Omega2000 and OmegaPrime. Data 
reduction was carried out using standard procedures in 
IRAF. Photometric calibration is based on the 2MASS 
point source catalog. 

3. RESULTS 

Fig.|21displays contours of SCUBA 850 /jto emission to- 
wards ISOSS J18364-0221 overlayed onto a near-infrared 
color composite. The IRAS and ISOSS FIR point source 
is clearly resolved into two emission peaks in the submil- 

limeter. The FWHM is 17" x 10.5" (PA 10°) for the 

eastern core (SMM1), while the western core (SMM2) is 
more extended. The fluxes are -Fsso/jm = 2.11Jy and 
^Wm = 13.5 Jy for SMM1 and F 850flm = 2.85 Jy and 
-F450Atm = 15.8 Jy for SMM2. In addition to the pho- 
tometric uncertainties of the submm fluxes, the 850 /im 
flux may be contaminated by the CO (3-2) transition. We 
checked the CO(3-2) contribution in the SCUBA 850 /ira 
band using our HHT CO (3-2) map and found the level to 
be within our overall calibrational uncertainty. Deriva- 
tion of the dust temperature from submm data alone 
heavily depends on the chosen dust emissivity index (3 
for which values between 1 and 2 are reported in the 
literature. Using an emissivity index of /3 = 2 we find 
Id = 14 K for SMM1 and T d = 12 K for SMM2. For lower 
values of (3 the temperatures would be higher. In order 
to constrain both (3 and the dust temperature we use the 
measured 100 /im flux (Fioo/jm = 57.7 Jy) of the whole 
ISOSS J18364-0221 cloud that places an upper limit for 
the emission of the individual cores. Using HIRES pro- 
cessed maps of the region, we find that the relative flux 
ratio of SMM1 to SMM2 at 100 ^m is approximately 10:1. 
Using this ratio we derive (3 = 1.6 and T d = 21 K for 
SMM1 and (3 = 1.8 and T d = 13.8 K for SMM2. Since 
this is an upper limit we take the average values for the 
rest of the paper: (5 = 1.8 ± 0.2, T d = 16.5^ K and 

(3 = 1.9 ± 0.1, T d = 12.8i L -gK for SMM1 and SMM2, 
respectively. 

We estimate the masse s of the cores using the rela- 
tion of lHildebranl lfl98l M d = F x d 2 /K X Bx{T) with a 
gas-to-dust ratio of 100. The mass of SMM1 is M = 
To ± 30 M Q , SMM2 has a mass of M = 28Ot^M . 
We used K 850Aim = 0.018 for SMM1 (gas+dust, thin 
icemantles, n ~ 10 6 cm~ 3 ) and for SMM2 Ksso^m = 
0.010 (gas+dust, thin i ceman tles) taking the values from 
lOssenkopf fc Hennind ilL99l L 

The submillimeter emission is highly correlated with 



the NIR extinction, as seen by the reddening and sur- 
face density of background stars in Fig. We have 
constructed a quantitative extinction map from our J, 

H, and Ks images accor ding to the Nicer method of 
lLombardi fc Alvesl Ipffll (Fig.EJ. The map reveals that 
the compact far-infared and submm sources are embed- 
ded in the center of a much more extended cloud complex 
which is also visible as an optical dark cloud Lynds 541. 

In contrast to the submillimeter maps which resolve 
individual features on typical scales of individual cloud 
cores (0.1 pc), the ISOSS and IRAS beams constrain the 
average temperature and mass of the whole dense cloud 
center on approximately 1 pc scale. The small chopping 
throw of SCUBA (2') results in spatial filtering of struc- 
tures that are comparable or larger in size. We account 
for an average column density of Ay ~ 7 (see Fig. to 
correct the integrated flux of the chopped SCUBA mea- 
surements. We fit a single modified black body to the 
data longwards of 100 /im (Fig. 0J) and derive a temper- 
ature of 14.6^'q K and (3 = 2. The total luminosity of 
ISOSS J18364-0221 is L ~ 800 L Q . 

The mass of the region is M = 900^3^ M© (^170^111 = 
25cm 2 g _1 ) with the uncertainty dominated by the de- 
rived average dust temperature. A lower limit for the 
mass can be computed from our extinction map, us- 
ing the relation M = (AOd) 2 /j^-^ j4y(i) given by 
IDickmanl l|1978|) . With a dust-to-gas ratio of A^ /iy = 

I. 87 x 10 21 cm" 2 mag -1 ijSavage fc Mathislll979f) we ob- 
tain a mass of M > 460 M for ISOSS J18364-0221 and 
M ~ 3200 Mq for the whole complex shown in Fig. EJ 

Ammonia (NH3) is an ideal tem perature probe for 
the molecular gas in dense r egions ijHariu et alJ 119931 
Wal mslev fc Ungerechtsiri983(l . From the relative inten- 
sity of the (J,K) = (1,1) and (2,2) inversion transitions of 
our ammonia (NH3) data, we derive a gas kinetic temper- 
ature of T kin = 11.6 ± 1.5 K towards SMM2 (Fig. [5). The 
inversion transitions (J,K) = (3,3) and (4,4) were not de- 
tected in neither of the cores (ctrms ~ 30 mK). The latter 
transitions have been used as sensitive tracers of HMCs 
and argue against the presence of such an object in ISOSS 
J18364-0221. SMM2 has a shallow density profile and is 
larger than the beam size of the ammonia measurements 
(40"). Therefore beam dillution is neglectable and we 
can directly compare the line and continuum measure- 
ments. The good match of Tki n and T d (in particular for 
(3 = 2) suggests thermal coupling between dust and gas, 
as it is expected for cold and dense cores. 

4. DISCUSSION 

We now discuss the nature of the two massive cold 
cores in ISOSS J18364-0221 and the implications for the 
initial conditions of high-mass star formation. 
SMM1 is cold (T w 16.5 K) and compact (effective ra- 
dius R ~ 0.2 pc). From our 450 //m data we derive a col- 
umn density of A^(H 2 ) = 2.2 x 10 23 cm -2 and a minimum 
central density of rt(H 2 ) = 1.4 x 10 6 cm~ 3 . Our measure- 
ment of optically thin H 13 CO + (3-2) transition towards 
SMM1 (Fig. E| provides an independent check for the 
temperature and (column) density of the high-density 
gas. Using the density und temperature from the dust 
continuum and appl ying the radiative transfer code from 
Schoi eret alJ l)2005f) . we find an abundance of [HCO/H 2 ] 
of 8 x 10 -10 which is in good agreement with the value 
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predicted for protostellar cores (Be rgin &: Langeilll997j) . 
Our molecular line measurements further indicate that 
the core is collapsing and drives an outflow: The optically 
thick H 12 CO + (3-2) transition shows a red-shifted self ab- 
sorption with respect to the optically thin H 13 CO + (3-2) 
line (Fig. EJ, sugg esting larg e scale infall towards the 
center of the core ilChoill2002t lEvanllT999T) . The CO(2- 
1) spectra show significant line wings around the two 
submm cores when compared with C 18 0(2-l), indicat- 
ing the presence outflows (see Fig. We derived the 
outflow para meters using ou r CO (2 -1) and C 18 0(2-l) 
data following iHenning et alJ l)2000bD . The mass of the 
outflow is estimated to Mf « 18 M Q and the dynami- 
cal timescale ta is ~ 1.8 X 10 4 years. The correspond- 
ing mass outflow rate is M « 1 x lO _3 M0/yr and the 
momentum supply rate F w 8.5 x 10 _3 Mq kms _1 /y r ) 
without correcting for the unknown inclination i of the 
outflow. This correction would further increase the de- 
rived values, which are already comparable to those ob- 
served for outflows in high-mass star f orming regions (e.g. 
iBeuther et al.l2002tlZhang et al.l2005[) . It must be noted 
however that both the red and blue lobes have complex 
morpholgies and higher resolution observations are re- 
quired to decide on the actual number of outflows driven 
by SMM1. No infrared counterpart is detected in our 
deep near-infrared images nor in the MSX mid-infrared 
data, which is however not surprising considering the ex- 
tinction of Ay > 120. With M ~ 75 M SMM1 would 
be massive enough to form an O type star and an associ- 
ated cluster. We therefore consider SMM1 as a promising 
candidate for a massive protostar during its initial gravi- 
ational collapse, i.e. in an evolutionary state comparable 
to the Class phase in low-mass star formation. 

SMM2 has very low gas and dust temperatures 
(T 12 K), moderate densities (central density n{B-2) ~ 
2.3 x 10 5 cm~3), and is quite extended (effective ra- 
dius R ~ 0.5 pc). Such conditions have been sug- 
geste d to characterize the earliest stage of star forma- 
tion ijEvans et alJ l2002). but have so far only been ob- 
served in low-mass prestellar cores. We investigate the 
dynamical state of the core considering the virial the- 
orem, with E mag + Spot = 2(E kin - E cxt ) being the 
condition for a gravitationally bound core. E mag de- 
notes the magnetic energy, E pot is the potential en- 
ergy, .Ekin the total kinetic energy, and £" e xt accounts 
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for external pressure. The kinetic energy can be writ- 
ten as E kin = E thcrm + E tulh = iNkT + \Mo\ urh , with 
cr 2 = AV 2 /(8 In 2) - kT/m, where AV = 0.90 km/s is the 
linewidth taken from our ammonia data (Fig. and m 
the mass of the NH3 molecule. The potentia l energy of 
a hom ogenous ellipsoid is E pot = 3GM 2 /5i? llLiliestroml 
IT99I . We find Epot = 9.5 x 10 38 J, E thcrm = 3.1 x 10 37 J, 
and E turD = 1.2 x 10 38 J. Neglecting the external pres- 
sure and the magnetic field energy the total kinetic en- 
ergy sums up to ~ 16% of the potential energy. We 
conclude that SMM2 is Jeans unstable. The small line 
width of the optically thin ammonia transitions sug- 
gests that the level of turbulence is not as high as in 
more evolved high-mass cores. There are two sources 
with very red colors towards SMM2 (H-K=3.7 and H- 
K=4.5) that are presumably low-mass Class I objects, 
based on their NIR and MIR fluxes with a spectral in- 
dex of a = dlog(AEx)/dlog(A) ~ 3. This classification 
is supported by the molecular outflow traced in CO(2-l) 
emerging from the young stellar objects (Fig. [2J. This 
supports the large-scale Jeans instability of the core and 
indicates that fragmentation and cluster formation has 
already started. Therefore SMM2 is a very promising 
candidate for a massive pre-stellar core, characterizing 
the earliest phase of st ar formation. 

Several authors fe.g. lEvans et a l. 2002) have suggested 
that the conditions during the earliest phases of high- 
mass star-formation should be similar to that prevailing 
in the regime of low-mass star birth. Our findings provide 
direct observational evidence that the initial conditions 
of high-mass star-formation are indeed characterized by 
very low temperatures and low levels of turbulence. 
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Fig. 1. — a) Incomplete 170 fim map constructed from slew-measurement by the ISOPHOT Serendipty Survey. ISOSS J18364-0221 
the compact source in the center, b) 100 fim map of the same field, produced by HIRES processing of the IRAS survey. 
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Right ascension 

Fig. 3. — Map of the visual extinction Ay towards ISOSS J18364-0221. Contours start at Ay = 4 and increase in steps of 4. The 
local ex tinction was determ ined from the reference region R. Colors for zero extinction were established using a simulated besancon model 
catalog I Robin ct al. 2003) towards the region. 
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Fig. 4. — Spectral energy distribution of ISOSS J18364-0221. For wavelengths longwards of 100 fim the SED is dominated by the optically 
thin thermal emission of large grains and well fitted with a modified black body with an emissivity index of /3 = 2, the point at 60 (im is 
not taken into account. 




